Transition to turbulence in particulate pipe flow 
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We investigate experimentally the influence of suspended particles on the transition to turbulence. 
The particles are monodisperse and neutrally-buoyant with the liquid. The role of the particles on 
the transition depends both upon the pipe to particle diameter ratios and the concentration. For 
large pipe-to-particle diameter ratios the transition is delayed while it is lowered for small ratios. A 
scaling is proposed to collapse the departure from the critical Reynolds number for pure fluid as a 
function of concentration into a single master curve. 

PACS numbers: 47.20.-k, 83.80.Hj, 83.50.Ha 
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More than a century after Reynolds' work pj, un- 
derstanding how turbulent regions grow in a pipe and 
bring the laminar Poiseuillc flow to fully developed tur- 
bulence is still not completely achieved. Above a crit- 
ical Reynolds number, the laminar flow is observed to 
be unstable, turbulent regions grow and are convected in 
the pipe. This flow regime is called intermittent. When 
the flow rate is further increased, the flow becomes fully 
turbulent |2j. In fact, the transition happens to be sub- 
critical and the flow is linearly stable for all flow rates 
0- A finite amplitude perturbation is needed to trigger 
the transition and the critical Reynolds depends upon 
its amplitude. For small perturbations, laminar motion 
is observed as far as Re ~ 10 5 but the transition in pure 
fluid can be reached for Re « 2100 provided that the 
perturbation is strong enough to allow the growth of tur- 
bulent "puffs" 0. Recent studies have investigated with 
different kinds of perturbations the nature of the unstable 
modes, either in the inlet region or in the fully developed 
flow [li]. 

The objective of the present work is to examine how 
transition to turbulence is affected by the presence of sus- 
pended particles in the simplest case of neutral buoyancy. 
More specifically, we focus upon determining the transi- 
tion threshold between the laminar and the intermittent 
regime as a function of the particle volume fraction <j> of 
the suspension. Because the particles are neutrally buoy- 
ant and largely drag-free, the present study is related to 
recent work which has examined global subcritical stabil- 
ity behavior of plane Couette flow forced by the presence 
of a single spherical bead or a spanwise wire 0, Q . This 
work also has a practical aspect as it is related to pipeline 
flow of slurries. 

Experiments are performed with four sets of 
spherical polystyrene particles having density 
p = 1.0510 ± 0.001 g.cm -3 and diameters d 
presented in table [fl To obtain neutral buoyancy, the 
densities of the fluid and of the particles are matched. 
We choose as a fluid a mixture of 22 % glycerol and 
78 % water by mass. The temperature of the mixture 
is maintained at 25 ± 1° C by using a thermostated 
bath as a fluid reservoir in the fluid circulating loop. 
At this temperature, the viscosity of the mixture is 
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TABLE I: Particle diameters and pipe to particle diame- 
ter ratios. The smaller particles were supplied by Kodak 
(Rochester, NY USA) while the others were from Maxiblast 
(South Bend, IN USA). 



H = 1.64 ± 0.03 cP. 

The experimental set-up consists of a straight and hor- 
izontal cylindrical glass tube of 2.6 m length mounted on 
a rigid support structure. Two different tubes having dif- 
ferent inner diameters D± = 8 mm and D 2 = 14 mm are 
used in the experiments. These tubes are longer than the 
entry lengths necessary for the laminar flow to fully de- 
velop at Re w 2000, L e (D{) ~ 0.6 m and L e (D 2 ) ~ 1 m. 
The pipe to particle diameter ratios, D/d, used in the ex- 
periments by combining the different particles and tubes 
are indicated in tabled In order to ensure that the flow in 
the pipe is undisturbed by perturbations from a pump, 
the flow is driven by gravity. The suspension is deliv- 
ered to the tube by overflow from a tank positioned at 
a fixed height to an outlet of variable height, passing 
through the glass tube. A Moineau progressing cavity 
pump (PCM model MR2.6H24) carries the suspension 
back to the overflowing tank, but is isolated from the 
flow through the glass tube. 

With these flow conditions, the transition for pure fluid 
was found to take place at Re w 2300, but the flow was 
very sensitive to any kind of perturbation applied to the 
pipe. In order to control the transition with a known 
perturbation, a ring (solid annulus) fitting tightly into 
the tube and of thickness 1.5 mm was inserted in the 
pipe at its entrance. The perturbation produced by the 
ring then controlled the transition and lowered it down 
to Re 2100, which is the lowest value that could be 
reached with this kind of perturbation. The influence of 
the particles on the transition has been studied with this 
ring. 

The different flow regimes, i.e. laminar, intermittent 
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or turbulent, are identified by visualization as well as by 
measuring the pressure drop between the entrance and 
the exit of the glass tube. 

A classical way to detect the pipe-flow transition for 
pure fluid is to inject dye and to observe the wavering of 
the streak as it passes down the pipe pIQ In the present 
experiment with particles, direct visualization of the par- 
ticle motion in the suspension is a very straightforward 
way of observing the birth and convection of "turbulent 
puffs" in the flow. Whereas motion in the laminar regime 
is characterized by parallel particle trajectories, motion 
in a turbulent puff presents a strong mixing in the radial 
direction. When the flow rate is increased, we observe 
more and more turbulent puffs in the flow, and finally 
these merge once the turbulent regime is established. 

To provide a more quantitative indicator of transition, 
we also measured the pressure drop between the entrance 
and the exit of the glass tube with electronic manome- 
ters (Newport Omega PX 154). Pressure fluctuations 
due to the larger pressure drop caused by the turbulent 
puffs could then be clearly identified with a signal ana- 
lyzer (Hewlett Packard 3562 A). Figure ^ shows pressure 
drop spectra in the low-frequency range for the laminar 
flow and at the onset of intcrmittency, both for a particle 
volume fraction <f> — 0.1. The spectrum in the intermit- 
tent regime displays strong fluctuations for frequencies 
under 0.1 Hz. Using the spectral signature, the onset 
of intermittency can thus be clearly identified even for 
large volume fractions where the visualization is not pos- 
sible. The spectrum of pressure fluctuations in the flow 
driven by the pump was observed to consist of "contin- 
uous noise" of the type shown by figure ^ phis a strong 
peak at the frequency of the pump rotation with its as- 
sociated harmonics. By contrast, the fluctuations when 
the flow is driven by gravity are observed to have only 
the "continuous noise". Furthermore, this noise was de- 
termined to be unchanged when the pump was off and 
the overflow tank was filled manually. Consequently, the 
spectrum observed in the gravity-driven flow is consid- 
ered to be independent of the pump. 

When the threshold of intermittency was reached, we 
measured the critical flow rate Q c and deduced the crit- 
ical Reynolds number Re c — AQcp/irfiD with the pipe 
diameter D as the length scale. To measure Q, we col- 
lected a given volume of the suspension without altering 
the flow, by capturing fluid at the outlet to the thermo- 
statically controlled reservoir. The time needed to ob- 
tain this volume yielded Q c . The particles collected in 
this volume were then sieved, rinsed, dried, and weighed 
to provide the mean particle volume fraction (f> in the 
suspension flow. Measurements of the critical Reynolds 
number associated with the start of intermittency have 
been carried out for different combinations of particles 
and tubes presented in table [I] and for different concen- 
trations <p ranging from to 0.3. 

Figure[2]displays the critical Reynolds number Re c as a 
function of the particle volume fraction <p for the different 
combinations of particles and tubes. The results indicate 
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FIG. 1: Fluctuations of the pressure drop across the pipe for 
the laminar and the intermittent regime in a suspension of 
particle volume fraction <f> = 0.1. The zero frequency value of 
the spectrum, related to the mean value of the pressure drop, 
is much larger and not represented on the graph 
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FIG. 2: Critical Reynolds number Re c as a function of the 
volume fraction <f) of the suspension for different combinations 
of particles and tubes. Tube Di with d = 215 /jm (o), 510 fim 
(V), 780/im (□). Tube D 2 with d = 40 ^m (x), 215 /xm (•), 
510 yum (T), 780 ^m (■). 



two different situations: 

(i) For particles with D/d > 65 (x and •), the tran- 
sition is shifted to larger Reynolds numbers. The 
critical Reynolds number Re c is a monotonically 
increasing function of <j>. Furthermore, the data for 
the two sizes collapse onto a single curve, so that 
at a given in that range of D/d, the transition 
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threshold does not depend on the diameters of the 
particles and of the tube. 

(ii) For particles with D/d < 65, the behavior depends 
on <fi. For small </>, the transition is moved to lower 
Reynolds number and Re c decreases with increas- 
ing (f>. The magnitude of the decrease and in par- 
ticular the minimum Re c reached as </> is increased 
depends on the diameters of the particles and the 
tube. In particular, for the 780 /im particles in 
pipe D\ for which D/d ~ 10 (□), intermittency 
is observed at Re c ~ 1000 for <j> ~ 0.07. This is 
the smallest Re c found in this work. For larger </>, 
Re c increases with increasing <j> and the transition 
is eventually delayed. It is worth noticing that the 
data for combinations of particles and tube hav- 
ing roughly the same D/d collapse onto the same 
curve (o and T for D/d ~ 27 — 37 and V and ■ for 
D/d~ 16- 18). 

The behavior for D/d > 65 seems simpler since it is 
independent of the diameters of the particles and of the 
tube. The observed delay to the transition is expected 
to be due to an enhancement of the viscosity caused by 
the suspended particles. Common models for the effec- 
tive viscosity /i e of a suspension do not depend on the 
size of the particles and give an increase of the viscosity 
with increasing average concentration <fi. A well-known 
example is Krieger's viscosity H,Ei, which expresses the 
effective viscosity y, e as a function of <fi according to the 
law: 
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where fi is the viscosity of the pure fluid and <f m = 0.68 
the random close packing (i.e. maximum) concentration 
for spherical particles. This empirical formula initially 
determined at low Reynolds number is commonly used 
at finite Reynolds. 

The critical Reynolds number of the suspension using 
equation ifTJl. Re cs = Re c [i/fi e , is plotted as a function 
of (j> in figure We observe that for D/d > 65 (x and 
•) Re cs is approximately independent of <j> and remains 
close to the value for the pure fluid, Re cs ~ 2100. For 
4> > 0.25, we observe a steep increase in this thresh- 
old (• and also ▼), suggesting the presence of an ad- 
ditional mechanism for dissipation beyond the viscosity 
enhancement observed in Stokes flow suspensions. How- 
ever, we have only a limited range of concentration data, 
because the large flow rates required to achieve transi- 
tion at the elevated effective viscosity result in a large 
pressure drop. The large pressure exceeds the range of 
the pressure gauges and also necessitates the direct use 
of the pump instead of gravity-driven flow that we have 
chosen to consider only in this study. 

This simple scaling with Krieger's viscosity is not suf- 
ficient to obtain a collapse of the curves for D/d < 65 
(see figure for o, V, □, T, ■). However, we can no- 
tice a similarity between these curves. At low <j>, there 
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FIG. 3: Critical Reynolds number of the suspension using 
Krieger's viscosity Re cs = i?e c /io//i e as a function of particle 
volume fraction <f>. Tube Di with d = 215 /im (o), 510 /tm 
(V), 780/xm (□). Tube D 2 with d = 40 ^m (x), 215 /xm (•), 
510 urn (▼), 780 ^m (■). 
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FIG. 4: Scaled difference between the critical Reynolds num- 
ber of the suspension using Krieger's viscosity Recs and that 
of the pure fluid Re c o — 2100, i.e. (Re cs — Re c o)D/d as a 
function of (j>D/d for <j> < 0.25. The values of the parame- 
ters are as in figure El (only the data for particles which cause 
earlier transition have been kept). 



is a first regime where Re cs decreases steeply with in- 
creasing <f> and the different curves appear to collapse. It 
should be mentioned that a very small decrease is also ob- 
served for the 215 /mi particles in pipe D2 (•) but it can 
be considered to be within error bars. Above a critical 
volume fraction <f> c which depends on the pipe to parti- 
cle diameter ratio D/d, there is a second regime where 
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the curves eventually reach minimum values of Re cs and 
remain approximately independent of <j) for larger vol- 
ume fractions. These minimum values of Re cs decrease 
with increasing D/d. Moreover, the data for combina- 
tions of particles and tube having similar values of D/d 
seem to collapse onto the same curve. These last obser- 
vations lead us to scale the difference between the critical 
Reynolds number of the suspension Re cs and that of the 
pure fluid Re c o — 2100 as well as 4> with D/d as dis- 
played in figure 0] This new scaling provides a collapse 
of the different curves for all D/d. The scaled difference 
(Re cs — Re c o)D/d initially decreases rapidly and approxi- 
mately linearly with <pD /d, then saturates and eventually 
increases slightly at larger <f>D/d. The results on figure 0] 
suggest <p c ~ d/D. 

These data were obtained with the ring at the entrance 
of the pipe but similar results were also observed without 
it. This suggests that the subcritical transition is trig- 
gered by the particles. A particle introduces fmctuational 
velocities whose form and coupling to the mean flow vary 
with the particle Reynolds number Re p — Re d 2 / D 2 . Un- 
like the case of uniform flow, the influence of Re p on the 
structure of the dual wakes caused by shear flow around 
suspended bodies is, while under study 0, Il4j , not well 
understood. It is plausible that the influence of particle 
size in reducing Re cs at a fixed <f> results from the increase 
in the particle scale Reynolds number with d/D. With 
increasing Re p , the disturbance flow caused by the par- 
ticle is presumably less efficiently dissipated by viscous 
action, thus allowing for stronger coupling to the bulk 
flow. 

A more refined analysis should take into account the in- 
fluence of particles on the velocity and concentration pro- 
file in the pipe flow. In conditions of very low Reynolds 
number and high concentration, particles migrate to- 
wards the center of t he p ipe and blunt the velocity pro- 
file, see for instance [l(| ■ This effect which would tend 
to reduce the effective viscosity in the pipe is suspected 
to be present for the 780 fim particles in pipe D\ for 
4> > 0.25 (see □ in figure [3J|. There exists a second type 
of migration, the so-called "tubular pinch effect" , which 
is inertial and causes a single particle to move to a po- 



sition at a distance of 0.31? from the axis [l2L IT3T |. This 
effect is observed with the larger particles in both pipes. 
The present study seems to suggest that the particles 
alter the threshold of the subcritical transition through 
coupling of the base flow to velocity fluctuations rather 
than the base flow itself through their migration but this 
requires confirmation. 

This work leads to definite conclusions regarding the 
influence of suspended neutrally-buoyant solids upon the 
transition away from laminar flow. The influence de- 
pends both upon the pipe to particle diameter ratios and 
concentration. For D/d > 65, neutrally-buoyant parti- 
cles cause a delay in transition to larger Re c . This effect 
can be explained by the enhancement of the effective vis- 
cosity of the suspension. However, for <f> > 0.25, the delay 
in transition is found to be substantially larger than can 
be explained by a simple renormalization by effective vis- 
cosities. For D/d < 65, the behavior is quite different and 
neutrally-buoyant particles alter the transition to turbu- 
lence in pipe flow to smaller Re c . Scaling the departure 
from the critical Reynolds number for pure fluid as well 
as the concentration with D/d gives a master curve for 
the transition for all D/d. The most plausible explana- 
tion for the reduction of the critical Reynolds number, 
though one still in need of confirmation, is that the fluc- 
tuations induced by the particles trigger the subcritical 
transition. 
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